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ABSTRACT 



Context. Bars are very common in the centre of the disc galaxies, and they drive the evolution of their structure. The state-of-the-art 
imaging and redshift surveys of galaxies allow us to study the relationships between the properties of the bars and those of their hosts 
in statistically significant samples. 

Aims. A volume-limited sample of 2106 disc galaxies was studied to derive the bar fraction, length, and strength as a function of the 
morphology, size, local galaxy density, light concentration, and colour of the host galaxy. The sample galaxies were selected to not 
be strongly disturbed/interacting. 

Methods. The bar and galaxy properties were obtained by analysing the r-band images of the sample galaxies available in the Sloan 
Digital Sky Survey Data Release 5. 

Results. The bars were detected using the ellipse fitting method and Fourier analysis method. They were tested and calibrated with 
extensive simulations on artificial images. The ellipse fitting method was found to be more efficient in detecting bars in spiral galaxies. 
The fraction of barred galaxies turned out to be 45%. A bar was found in 29% of the lenticular galaxies, in 55% and 54% of the early- 
and late-type spirals, respectively. The bar length (normalised by the galaxy size) of late-type spirals is shorter than in early-type or 
lenticular ones. A correlation between the bar length and galaxy size was found with longer bars hosted by larger galaxies. The bars 
of the lenticular galaxies are weaker than those in spirals. Moreover, the unimodal distribution of the bar strength found for all the 
galaxy types argues against a quick transition between the barred and unbarred statues. There is no difference between the local galaxy 
density of barred and unbarred galaxies. Besides, neither the length nor strength of the bars are correlated with the local density of 
the galaxy neighbourhoods. In contrast, a statistical significant difference between the central light concentration and colour of barred 
and unbarred galaxies was found. Bars are mostly located in less concentrated and bluer galaxies. 

Conclusions. These results indicate that the properties of bars are strongly related to those of their host galaxies, but do not depend 
on the local environment. 

Key words, galaxies: evolution — galaxies: fundamental parameters — galaxies: spiral — galaxies: structure — galaxies: statistics 



1. Introduction 

Strong bars are observed in optical images of roughly half of all 
the ne a rby disc galax ie s dMarinova & JogeeL l2007t iReese et all 
l2007t IBarazza et all I2008I) . This fraction ri ses to about 
70% when near-infrared images are analysed (|Knapen et ail 
2000; lEskridge et all l2000t iMenendez-Delmestre et all I2007I) . 
The presence of a bar can be found by visual inspection 
(e.g., de Vaucouleurs et ail Il99ll hereafter RC3), by analysing 
th e shape an d orien t ation of the galaxy isophotes (e.g._ 



Bara zza et alll2008l). although there are some claims that this is 
not the case dShefh et aTll2008l) . 



Wozniak et al 



. .19951; [Laine et all 120021; iMarinova & JogeeL 
2001: I Barazza et al I 120081) . or by stu dying the Fourier 
modes of the light d i stribu ti on (e.g.. lOhtaetall \l 990 ; 
Elmegreen & Elmeg reenl 119851 ; [Aguerri etall 119981 12000a). 



Therefore, bars are a common feature in the central regions 
of disc galaxies of the local universe. But, this is also true at 
high redshift. In fac t, the bar fraction apparently remains con- 
stant out to z « 1 (Jo g ee et all |2004; El megreen et all [2004; 
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The presence of bars in the centre of lenticular and spiral 
galaxies make them ideal probes of the dynamics of the cen- 
tral regions of discs. In fact, bars are efficient agents of angu- 
lar momentum, energy, and mass redistributio n. They act on 
both luminous and dark mat t er components dWeinbergj 1 1985; 
Debattista & Sellwoo ! 119981 12000t [Athanassoulal 120031) driv- 
ing the evolution of galaxy structure and morphology. In par- 
ticular, the amount of angular momentum exchanged is re- 
lated to specific properties of the galaxies, su ch as the bar 
mass, halo density, and halo velocity dispersio n ( Athanassoula, 
120031; ISellwoodl 120061; ISellwood & Debattistal 120061) Moreover, 
they funnel material t owards the galaxy centre buildi ng bulge- 
like structures (e.g., Kormendy & Kennicutt, 2004), nuclear 
star-f orming rings (e.g., iButa et ail [2003), and nuclear bars 
(e.g., Erwinl 20041). a nd feeding the central black hole (e.g. 



Shlos man et all l2000t) 
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Bars play an important role in bulge formation. Major 
mergers or monolithic collapse are the c lassical theo ri es for 
bulge formation in disc ga laxies (e.g., Eg gen et ail, 119621; 
Kormendy & Kennicutt, 2004). But bulges can be built via mi- 



1990; 



2004 



nor mergers (I Aguerri et all 12001 at lEliche-Moral et all [2006) 
or secular evolution processes produced by bars. In fact 
according with the results of N-body simulations, the in- 
ner parts of a bar inflate after a few bar rotations be- 
cause of large-scale violent bending instabilities and set- 
tle with an increased thickness and vertical veloci ty dis- 
persion (e. g., Combes & Sanders l Il98lt IC ombes et al 
Rah a et all 1 199 it lAthanassoulal 120031; IDebattista et al 
Athanassoulal l2005t iMartinez-Valpuesta et all 120061) . This 
leads to the establishment of the connection between the 
bar-buckl ing mechanism and th e form at ion of boxy /peanut 

bulge s ([Bureau & Athanassoula, Il999t iBureau & Freeman! 

1999; IChung & Bureaut 120041) The buckling instability does 
not destroy the bar and forms a central stell ar condensation 
remin i scent of the bu l ges of late-type spirals ( IDebattista et all 
2004; Athanassoula, 2005), in agreement with early find- 
ings by iHohll (Il97ll) . Observational evidences of secu- 
lar bulge formation includes the near-exponential surface 



brightness pro fi les of some bu l ges ( Andredakis et all 



1996 : Ide Jongl Il996t jCarollo et all 



1995 



2001 



2005 



Courteau et al , 

Pricto ct all 1200 It iMacArthur ct all 120031 lAguerri et all 
Mendez-Abreu et all l2008at iFisher & Drorvt 120081). a corre- 
lation between bulg e and disc scale length dMacArthur et all 
2003; lAguerri et all l2005t iMendez-Abreu et all l2008al). the 
similar colours of bulges and inner discs dPeletier & Balcellst 
1996: lCourteau et all [19961 ICarollo et alll2007l) . substantial ro- 
tation (iKormendv & Kennicuttl 120041) . and the p resence of B/P- 
shape d bulges in w 45% of edge-on galaxies dLiitticke et all 
2000). Recently, the connection between B/P-shaped bulges 
and bars has also been con firmed in face-on barred galaxies 
dMendez-Abreu et alll2008cT) . 

The most important parameters of bars are the length, 
strength, and pattern speed. Their evolution depends on the ef- 
fectiveness of the angular momentum exchange between lumi- 
nous and dark matter. Different methods have been proposed to 
measure bar properties. 

The bar length can be obtained by eye estimates on galaxy 
images (Kormendy, Il979t iMartint 11995), locating the maxi- 
mum ellipticity of the galaxy isopho t es dWozniak et ail j 1995; 
lLaine et ail [2002; Mari nova & Jogeel 120071). looking for varia- 
tions of the isophotal position angle ( Shet hetail 120031; lErwiri 
120051) or of the phase a ngle of the Fourier modes of the 



galax y light distribution (Ouill enet all 1 19941 lAguerri etai 



2003 ). analysing the bar-interbar contrast dAguerri et all 12000a , 
2003), or by photom etric decomposit ion of the surface 



brightness distribution dPrieto et all 119971; lAguerri et ail [2005; 
Laurikain en et all 12005 ). The previous techniques reported that 
the typical bar length is about 3-4 kpc, and is correlated with 
the disc scale, sug gesting tha t the two components are affect- 
ing each other dAguerri et all [20051 IMarinova & Jogeel l2007t 



lLaurikainen et al. , 2007). 



The bar strength c an be derived by measu ring the bar torques 
dButa & Block! 1200 1|) . isophotal ell i pticity (|Martinet & Friedli[ 
1997t lAguerri 119991; IWhvte et all 120021 IMarinova & J ogee. 



2007T). the maximum amplitude of the m — 2 Fourier m ode 



(Atha nassoula & Misiriotist l2002t lLaurikainen et all 120051). or 
integrating the m — 2 Fourier mode in the bar region (Ohta et al., 
1990; Aguerri et al., 2000a). The bar strength is almost con- 
stant with Hubble type dMarinova & Jogeel 120071). but lenticula r 
galaxies host weaker bars than spirals (Laurikainen et aLll2007l) . 



The pattern speed of bars can be indirectly estimated by 
identi fying rings with th e locat i on of the Li n dblad resonances 
(e.g., IVega Beltran et all Il997t iJeong et all 120071) . matching 
the observed velocity and density fields with numerical mod- 
els of t he gas flows (e.g., Lindblad et all 19961 Aguerri et all 
l2001bl: IWeiner et all 1200 lb iRautiainen et all 120081) . analysing 
the offset and shape of the dust lanes, which trace the lo - 
cation of shocks in the gas flows (e.g.. lAthanassoulal 1 1992b . 
looking for colour changes (Aguerri et all t 2000a) and min- 
ima in the star formation dCepa & Beckmanl 1 19901) outside the 
bar region, or by adopting the Tremaine-Weinberg method (see 
ICorsinit l2008t for a review). Th e last is a model-independen t 
way to measure the pattern speed dTremaine & Weinbergtl 19841). 
which was successfully a pplied t o singl e (Merrifield & Kuiiken, 
119951 IDebattista et all l2002l lAguerri et all 120031). double 
(ICorsini et all 120031) and dwarf barred galaxies dCorsini et all 
120071) too. Observed pat tern speeds imply that barred galax - 
ies host maximal discs (IDebattista & Sellwoodl Il998l |2000), 
since bars in dense dark matter halos are rapidly deceler- 



ated by dynamical friction (Weinbergl 1 19851 ISellwoodl l2006t 
ISellwood & Debattistal 120061) . 

Nowadays, the large galaxy sur veys, such as the Sloan 
Digital Sky Survey d York et all feOOOl hereafter SDSS) allow us 
to study the bar properties in samples of thousands of galax- 
ies, obtaining for the first time statistically significant results. By 
studying a volume-limited sample of ~3000 galaxies in the local 
universe extracted from the SDSS, we plan to address three main 
issues: the first is to assess, by means of extensive tests on sim- 
ulated galaxies, the advantages and drawbacks of the two more 
common methods used detecting bars: the ellipse fitting and the 
Fourier analysis. The second is to investigate the possible differ- 
ences of the bar fraction and bar properties (length and strength) 
with the morphological type, ranging from SO to late-type spi- 
rals. The third is to understand how the properties of the host 
galaxies affects the formation of the bar and its properties. 

The paper is organised as follows. The galaxy sample is pre- 
sented in Sect. [2j the methods we adopted to detect the bars in 
the sample galaxies are explained in Sect.[3j they are tested using 
artificial galaxy images in Sect. [4j the fraction of bars are given 
in Sect. [5J the bars properties (length, and strength) are reported 
in Sect. [6] and compared to galaxy properties in Sect. [71 con- 
clusions are given in Sect. [8] Throughout this paper we assume 
// =100kms-' Mpc" 1 . 

2. Sample selection and data reduction 

The sample galaxies were selected in the spectroscopic 
catalogue of the SPSS Da ta Release 5 (SDSS-DR5, 
lAdelman-McCarthv etai! 120071) . From the -675,000 galaxies 
available in the catalogue, we took all the galaxies in the redshift 
range 0.01 < z < .04 and down to an absolute magnitude 
M r < -20 (~ M*, iBlanton etail [2005). This represents a 
volume-limited sample, because the apparent magnitude of a 
galaxy with M r = -20 at z = 0.04 (m r ~ 15.5) is within the com- 
pleteness limit (m r - 17.77) of the SDSS-DR5 spectroscopic 
catalogue. 

In order to deal with projection effects, we restricted our 
sample to galaxies with b/a > 0.5, a and b being the semi-major 
and semi-minor axis lengths of the galaxies. For disc galaxies, 
this is equivalent to say that we have selected objects with in- 
clination i < 60°. Although the cut in the observed axial ratio 
of the sample galaxies introduces a bias in the selection of the 
elliptical ones, this will not affect the results of the paper since 
we are interested only in the properties of disc galaxies. 
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Then, we rejected all galaxies with close neighbours. To this 
aim, we excluded all the galaxies with a companion which was 
closer than 2 x rgo, where rgo is defined as the radius which 
contains 90% of the total galaxy light. In addition, the compan- 
ion must be within + 3 magnitudes with respect to the magni- 
tude of the target galaxy to be excluded. In this way, galaxies 
with faint companions or possibly contaminated by faint fore- 
ground/background objects are not discarded in our study. The 
resulting sample co nsisted of 3 060 gala xies. 

According to Mari nova & Jog ee (2007) the bar length of 
barred galaxies in the local universe ranges between about 0.5 
and 5 kpc, with a mean value ~3.5 kpc. At z = 0.04 the shortest 
bar length (rb ar = 0.5 kpc) projects onto 2.17 pixel (0786) in the 
SDSS images, which have a scale of 073946 arcsec pixel -1 . The 
PSF of the SDSS images can be mod elled assuming a Moffat 
function (see e.g. iTruiillo et al.L 12001). We fit a bidimensional 
Moffat function of the form 



Table 1. Median values of the light concentration and central 
velocity dispersion for the different galaxy types. 



PSF(r) = 



KT- 



(i) 



to several stars in each galaxy field obtaining a typical FWHM 
and /? parameter of 1709 (2.77 pixel) and 3.05, respectively. This 
means that the bars with a length of 0.5 kpc are not resolved. 
According to the tests on artificial galaxies we performed (see 
Sect. 4), the smallest bars that we are able to recover in the SDSS 
images have a length of ~ 9 pixel. This corresponds to ~ 0.5 kpc 
at z = 0.01 and ~ 2 kpc at z = 0.04. Therefore, a value of 2 
kpc is a more reliable limit on the actual resolution of the bar 
length throughout the range of distances covered by our sample 
galaxies. 

The visual morphological classification given by the RC3 
was available only for a subsample of 612 galaxies. Automatic 
galaxy morphological classifications divide galaxies according 
to some photometric observables. In particular, the light con- 
centration is s t rongly correlated with the Hubble type (e.g . , 
lAbraham et all 119961: IConselice et all 120001: IConselicel [2003). 
In fact, it is greater in early- than in late-type galaxies. We de- 
fined the light concentration as C — rgo/rso, where r$o and 
are the radii enclosing 50% and 90% of the total galaxy light, 
respectively. These radii are available in the SDSS database for 
all objects of our sample. We calculated the median values of C 
for the ellipticals (T < -4), lenticulars (-3 < T < -1), early- 
type spirals (0 < T < 3), and late-type spirals (T > 4) of the 
RC3 subsample. They are given in Table [TJ and shown in Fig. [T] 
Although the dispersion of the data is large, the median value of 
the light concentration decreases with increasing Hubble type. 

If the light traces mass, then we expect a relation between 
the mass concentration and morphological type too. The cen- 
tral mass concentration of a galaxy can be traced by the central 
velocity dispersion <x . This was available in the SDSS for 298 
of the sample galaxies listed in RC3. We calculated the median 
values of cr for the different Hubble types, and they are given 
in Table [TJ and shown in Fig.Q] The median velocity dispersion 
decreases with increasing Hubble type, too. 

Fig EJshows the relation between C and cr for those galax- 
ies in our sample with both the photometric and kinematic in- 
formation. Each galaxy can be assigned to a morphological bin 
according to its values of C and cr . They were assigned to the 
bin with closest median values listed in Table [TJ by minimizing 
the following equation: 



Galaxy Type 


C 


0"o 






(kms- 1 ) 


E 


3.18+0.15 


242+40 


SO 


3.10+0.31 


209+43 


SO/a-Sb 


2.53±0.50 


149+34 


Sbc-Sm 


2.10±0.24 


126±27 



Table 2. Comparison between our and RC3 classification. 





SO (our) 


Early-type spirals (our) 


Late-type spirals (our) 


E (RC3) 


0.16 


0.01 


0.01 


SO (RC3) 


0.53 


0.05 


0.03 


SO/a-Sb (RC3) 


0.26 


0.67 


0.37 


Sbc-Sm (RC3) 


0.05 


0.27 


0.59 



where C and cr' Q are the median values of C and ctq reported in 
Table [TJ 

Galaxies, for which only the light concentration was avail- 
able, were assigned to the morphological bin corresponding to 
the closest median value of C. 

We found that 26%, 29%, 20%, and 25% of the selected 
galaxies turned to be ellipticals, lenticulars, early-type and late- 
type spirals, respectively. In this work, we focused on the 2166 
disc galaxies which include the lenticulars and the early-type 
and late-type spirals. This represents our final sample. Among 
the disc galaxies 39%, 28%, and 33% were classified as lenticu- 
lars, early-type and late-type spirals, respectively. Table[2]shows 
the comparison between the RC3 morphological classification 
and our automatic classification for the sample galaxies listed in 
RC3. It lists the fraction of elliptical, SO, SO/a-Sb and Sbc-Sm 
galaxies as classified by RC3 among the galaxies we classified 
as lenticular, early- and late-type spirals. It is worth noticing that 
for all the disc galaxies more than 50% of the objects are as- 
signed to the same morphological bin by both the RC3 and our 
automatic classification. The better agreement is shown by the 
early-type spiral galaxies. Nevertheless, the contamination be- 
tween the different classes could be strong. 

The /--band image of each galaxy was retrieved from the 
SDSS archive. All the images were bias subtracted, flat-field cor- 
rected, and sky subtracted according to the associated calibration 
information stored in the Data Archive Server (DAS). 

3. Methods for detecting and analysing bars 

Three main methods have been proposed for detecting bars 
and analysing their propert i es (see for a review lErwinLl2005l: 
Michel-Dansac & WozniaH 120061) . They are based on the 
ellipse fit o f the galaxy isoph o tes dW ozniaket al., 1995 
Kna pen et all l2000t [Laine et all 120021: ISheth et all 12003 



Elmegreen et all 12004 : Jogee et all 20041 : Marinova & Jogee , 
2007; iBarazza et al .1 120081) . Fouri er analysis of the az- 
imuthal lumino s ity profile ([E lmegreen & Elmegreen, 

l et all 



1985; lOhtaetall Il990t lAguerrietall l2000at lButa 



d= yj(C - C') 2 + (logicro) - bgioQ) 2 , 



(2) 



2006J; lLaurikainen et al 1 120071) . and dec omposition of the 
galaxy surface-brightness distribution dPrieto et all [2001; 
lAguerri et ail l2005t jLaurikainen et aHl2005l) . 

In the present work, we have developed a fully automatic 
method for classifying barred and unbarred galaxies using the 
ellipse fitting (Sect. 13. It and Fourier analysis (Sect. 13.21 ). This 
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Fig. 1. Values of the light concentration C (top panel) and cen- 
tral velocity dispersion <x (bottom panel) as a function of the 
morphological parameter T. Only the sample galaxies in RC3 
are plotted. Median values of C and <x for ellipticals (diamond), 
lenticulars (asterisk), early-type spirals (triangle), and late-type 
spirals (square) are shown. 
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Fig. 2. Values of the light concentration C and central veloc- 
ity dispersion cr for the sample galaxies with cr available in 
SDSS-DR5. They are colour coded according to our morpholog- 
ical classification: the red, magenta, green, and blue dots corre- 
spond to the ellipticals, lenticulars, early-type and late-type spi- 
rals, respectively. Median values for the ellipticals (diamond), 
lenticulars (asterisk), early-type spirals (triangle), and late-type 
spirals (square) plotted in Fig. 1 are also shown. The full lines 
represent the locus of equal distance between different galaxy 
morphological types. 



automatic procedure has several advantages with respect to the 
user-dependent visual classifications. In fact, it is reproducible, 
it can be implemented and applied to large data sets. We first 
tested it by means of extensive simulations on a large set of arti- 
ficial disc galaxies (Sect. |4]l and then we applied it to the images 
of the sample galaxies. In a forthcoming paper we will obtain the 
structural parameters for the bulge, disc, and bar of all the sam- 
ple galaxies by applying a photometric decomposition of their 
su rface-brightness distrib ution based on the technique developed 
bv lMendez-Abreul d2008bl) . 



3. 1 . Ellipse fitting method 

This method is based on the fit of the galaxy isophotes by el- 
lipses. 

As a first step, each image was cleaned of field stars and 
galaxies. This was done by rotating the image by 180° with re- 
spect to the galaxy centre. Then, we subtract the rotated frame to 
the original one. The residual image was sigma-clipped to iden- 
tify all the pixels with a number of counts lower than 1 cr, where 
<x is the rms of the image background after sky subtraction and 
calculated in regions free of sources and far from the galaxy to 
avoid contamination. The value of the deviant pixels was set to 
zero. Finally, the clipped image was subtracted to the original 
one to get the cleaned and symmetrized image to be used in the 
analysis. 

The ellipses were fitted to the isophotes of the cleaned and 
symmetrized i mages of the 2 1 66 di sc galaxies using the IRAlfl 
task ELLIPSE dJedrzeiewskil 1 1987b . In order to get a good fit 
at all radii out to an intensity level corresponding to the back- 
ground rms , we i mple mented the fitting metho d described by 
iJogee et alj d2004l) and lMarinova & Joged d2007l) . This is an it- 
erative wrapped procedure, which runs the ellipse fitting several 
times changing the trial values at each fit iteration. At each fixed 
semi-major axis length, the coordinates of the centre of the fit- 
ting ellipse were kept fixed and corresponded to those of the 
galaxy centre. This was identified with the position of the inten- 
sity peak. The trial values for the ellipticity e and position angle 
PA were randomly chosen between and 1 and between -90° 
and 90°, respectively. The fitting procedure stopped when either 
the convergence was reached or after 100 iterations. The ellipse 
fit failed for a small fraction (60/2166) of the sample galaxies. 
For all the galaxies with properly fitted isophotes, we obtained 
the radial profiles of e and PA of the fitted ellipses. 

The ellipticity radial profile in a bright and inclined unbarred 
galaxy usually shows a global increase from low values in the 
centre to a constant value at large radii. At large radii the PA 
radial profile is constant too. The constant values of e and PA 
(on large radial scales) are related to the inclination and orien- 
tation of the line of nodes of the galactic disc. On the contrary, 
barred galaxies are characterised by the presence of a local maxi- 
mum in the elli pticity radial profile and constant PA in the b ar re- 
gion (see e.g JWozniak et al.Ul995tlAguerri et alll2000bl) . This 
is du e to the shape and orientation of the stellar orbits of the bar 
(see lContopoulos & Grosbollll989l : lAthanassoulalll992l) . 

This allowed us to identify bars by analysing the radial pro- 
files of e and PA. We considered that a galaxy hosts a bar when: 
(1) the ellipticity radial profile shows a significant increase fol- 
lowed by a significant decrease (Ae > 0.08), and (2) the PA 
of the fitted ellipses is roughly constant within the bar region 
(APA < 20° fl 

The values adopted for Ae and APA were determined 
by applying the method to artificial galaxies (see Sect, \4.2j . 
The bar length was derived as the radi us rj; at wh i ch the 
maximum e lliptic i ty was reached (e. g ., [W ozniak e t~afl 11995b 
lLaine et all |2002; Mari nova & Jogeel 120071) or as the radius 
r?£ at which the PA changes by 5° with resp ect to the value 
corresponding to the maximum ellipticity (e.g, IWozniak et all 



1 IRAF is distributed by NO AO, which is operated by AURA Inc., 
under contract with the National Science Foundation. 

2 We have considered that a global maximum in the ellipticity profile 
is produced by a bar when it is located at more than 3.5 pixels from the 
galaxy centre. This radial distance correspond to more than 3 times the 
FWHM/2 of the images. 
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19951 ISheth et all 12001 lErwin & Soarkel l200l lErwhl 12001 



Michel-Dansac & W ozniak. 2 0061) . 



This method has been already successfully applied by differ 
ent authors both in the optica l and near- infrared wavebands to 
detec t bars in galaxies at low ( Knapen et all l2QO0b lLaine et af 



2002| iMarinova & JogeeL 20071) and high redshift ( Jogee et al 
2004t lElmegree n et all 120041: ISheth et all 120081) . 



3.2. Fourier analysis method 

An alternative way to detect and characterise bars is with a 
Fouri e r analysis of the azimutha l lumin o sity prof ile ( Ohta et all 
1990; lElmegreen & Elmegreenl 119851: I Aguerri et aim OOOa). 
The method starts by deprojecting the image of each sample 
galaxy by a flux-conserving stretch along the minor axis by the 
factor 1 / cos i, where i is the galaxy inclination. Therefore, we 
needed a good estimation of the inclination and major-axis PA of 
the galaxy. They were derived from the ellipticity and PA of the 
ellipses fitted to the five outermost isophotes, where the total lu- 
minosity is dominated by the disc contribution. For each galaxy, 
the intensity of the outermost isophote corresponds to a value of 
lcr, where <x is the rms of the sky-subtracted background of the 
image. 

The deprojected luminosity profile, I(R, (/>), where (R, <f>) are 
the polar coordinates in the galaxy frame, is decomposed into a 
Fourier series 

I(R, 0) = -^-i + J] (A m (R) cos(m^) + BJR) sin(m0)) , (3) 



where the coefficients are defined by 
I(R, <p) cos(m<f))d(f), 



A,„(R) = 


n Jo 


and 




B m (R) = 


l f 

x Jo 



In 



I(R, <p) s,m(m(p)d<p. 



(4) 



(5) 



The Fourier amplitude of the m-th component is defined as 

A (R)/2 m = 0, 



UR) = 



y/A?„(R) + BUR) m * 0. 



(6) 



The even (m = 2,4,6, ...) relative Fourier amplitudes I m /Io 
of galaxies with bars are large, while the odd (m = 1,3,5,...) 
ones are small. In particular, the bar is evidenced by a strong 
m = 2 component. Similarly to the ellipse fitting method, we 
considered that a galaxy hosts a bar when: (1) the m = 2 relative 
Fourier component shows a local maximum (AfJz/Io) > 0.2), 
and (2) the phase angle of the m = 2 mode <p2 is roughly con- 
stant within the bar region (A02 < 20°). The values adopted for 
Mh/Io) and A02 were also determined by applying the method 
to artificial galaxies (see Sect. I4.31 >. 

The bar length r^°" n er was calc u lated u sing the bar/interbar 
intensity ratio as in [Aguerri et al.l (l2000 a). The bar intensity, 
4, is defined as the sum of the even Fourier components, 
Io + I2 + I a + while the inter-bar intensity, Ijb, is given by 
In - h + I a - h dOhtaetall Il990t | Elmegreen & Elmegreenl 
1985; lAguerri et all l2000al) . lohtaetalj d 19901) arbitrarily de- 
fined the bar length as t he oute r radius for which lb/lib = 2. 
However, Ag uerri et al.l (l2000al) pointed out that a fixed value 
of lb/ lib cannot account for the wide variety of bar luminosi- 
ties present in galaxies. Instead, they defined the bar length as 
the FWHM of the curve of Ib/Irmib - This method was applied 
bvlAthanassoula & Misi riofis] d2002l) to analytic models demon- 
strating its accuracy in measuring the bar length. 



4. Test on artificial galaxies 

4.1. Structural parameters of the artificial galaxies 

Extensive simulations on a large set of artificial disc galaxies 
were carried out to test the reliability and accuracy of the ellipse 
fitting (Sect. |3.1t and Fourier analysis (Sect. |3.2t in detecting 
bars. Moreover, they were used to fine tune the free parameters 
of the two methods, i.e., Ae and APA in the ellipse fitting, and 
Mh/Io) and Afc in the Fourier analysis. The surface-brightness 
distribution of the artificial galaxies was assumed to be the sum 
of the contributions o f three structural components: a bulge, a 
disc, and a bar (e.g. iPrieto et all I2001I; lAguerri et ail 120031 
120051; lLaurikainen et alT T2005). The surface-brightness distribu- 
tion of each individual component was assumed to follow a para- 
metric law, which has to be strictly considered as an empirical 

fitting function. 

The Sersic law dSersid fl968h was assumed for the radial 
surface-brightness profile of the bulge 



^bulgeW = ^0 



bulge 



10 



(7) 



where r e , /o.buige, and n are the effective (or half-light) radius, the 
central surface brightness, and the shape parameter describing 
the curvature of the profile, respectively. The value of b„ is cou- 
pled to n so that half of the total flux is a lways within r e an d can 
be approximated as b„ = 0.868/1 -0.142 dCaon et allll993l) . The 
total luminosity of the bulge is given by 



2 r(2n) 

^bulge = 27r/o,bulge(l ~ Qiulge)"^ - , 



(8) 



where ebuige is the observed ellipticity of the bulge and T is the 

Euler gamma function. 

The exponential law (Freeman! 1 1970b was assumed to de- 
scribe the radial surface-brightness profile of the disc 



Idisc(r) = /o,dis C e 



■r/h 



(9) 



where h and /o ; disc are the scale length and central surface bright- 
ness of the disc, respectively. The total luminosity of the disc is 
given by 



Ldisc - 2^/o ; disc(l - fdisc)^ , 



(10) 



where edisc is the observed ellipticity of the disc. 

Several parametric laws have been adopted in literature 
to describe the surface-b rightness dis tribution of bar s. Ferrers 
dLaurikainen et all 120051) . Freeman dFreemanl 1 19661) . and flat 
bars C Prieto et all 1 19971) were considered for the artificial galax- 
ies. 

The surface-brightness distribution was assumed to be 
axially symmetric with respect to a generalised ellipse 
dAtha nassoula et all 1 1990b . When the principal axes of the el- 
lipse are aligned with the coordinate axes, the radial coordinate 
is defined as 



y 



(1 " Cbar) 



1/c 



(ID 



where ebar is the ellipticity and c controls the shape of the 
isophotes. A bar with pure elliptical isophotes has c = 2. It is 
c > 2 if the isophotes are boxy, and c < 2 if they are discy. The 
parameters eb ar and c are assumed to be constant as a function of 
radius. 
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The radia l surface-brightness profile of a Ferrers ellipsoid 
dFerrersl,["l877l) is given by 



Mi-(£) 2 ) 





2\"bar+0.5 



r < r bar , 
r > r hM , 



where /o,bar, /"bar, and n bar are the central surface brightness, 
length, and a shape parameter of the bar, respectively. The to- 
tal luminosity is given by 



r Ferrers _ n„r ,A 



f 

JO 



(13) 



where e bal - is the ellipticity of the bar. The simu lated bars were 

[ enera ted by adopting n bal - = 2 following Laurik ainen et al.1 
2005). In this particular case, the total luminosity of the bar is 
given by 

I2T = *W1 - ^r)^ar^|- (14) 

The radial surface-brightness profile of a Freeman bar is 



bar 



/n,bar A/ 1 — ( 



>'ba 



(15) 



where /o.bar and r bar are the central surfa ce brightness and length 
of the bar, respectively dFreemani ri966). The total luminosity of 
a Freeman bar is 



t Freeman _ _/i \r „2 

Har ~ ^ e bal -)/o,bar^ bar , 



(16) 



where e bar is the ellipticity of the bar. 

Finally, the radial surface-brightness profile of a flat bar is 



rflat _ j 



1 



1 +, 



(17) 



where / ,bar and r bar are the central surface brightness and length 
of the bar. For radii larger than r b ar the surface-brigh tness pro- 
file falls off with a scale length r s dPrieto et all [1997). The total 
luminosity is 



LZ = -2 tt 7 ,bar (1 - ebar) r 2 s Li 2 (-e^ ,r -) , 



(18) 



where e bar is the ellipticity of the bar and Li 2 is the dilogarithm 
function (also know as the Jonquiere function). 

We generated a set of 8000 images of artificial galaxies with 
a Sersic bulge and an exponential disc. Among these galaxies, 
2000 have a Ferrers bar, 2000 a Freeman bar, 2000 a flat bar, and 
2000 do not posses a bar. 

The apparent magnitudes of the artificial galaxies were ran- 
domly chosen in the range 



10 < m r < 16, 



(19) 



corresponding to that of the sample galaxies. To redistribute the 
total galaxy luminosity among the three galaxy components, the 
bulge-to-total L\, u \ ge /L tot and bar-to-disc L bar /Ldisc luminosity ra- 
tio were taken into account. They were considered to be 



< /-/bulge //-tot < 0.7, 

and 

< Lbar/Z/disc < 0.3, 



(20) 



(21) 



following Laurikai nen et alj d2005l) . The adopted ranges for the 
effective radius of the bulge, scale-length of the disc, and bar 
length were sele cted according to the values measured for spi- 
ral ga la xies by iMollenhoff & Heidtl d2001l). iMacArthur et al.1 



(12) 120031). Mollenhoff (2004), LaurikainenetalJ (|2007|), and 



Mendez-Abre u et : alj d2008a[) . They are 
0.5 < r e < 3 kpc, 
1 < h < 6 kpc, 
and 

0.5 < r bar < 5 kpc, 



(22) 
(23) 

(24) 



respectively. The ellipticities of the structural components were 
also selected to mimic those measured in real galaxies (e.g, 
iMarinova & JogeeL 12007). They are 



0.8 < 1 - e bulge < 1, 
0.5 < 1 - e disc < 1, 
and 

0.2 < 1 - e bar < 0.7, 
with 

e b ulge ^ £d[ sc < 6 bar . 



(25) 
(26) 

(27) 
(28) 



Finally, the position angles of the three components were se- 
lected randomly between 0° and 180° to allow each component 
to be independently oriented with respect to the others. 

In each pixel of the resulting images noise was added to 
yield a signal-to-noise ratio (S IN) similar to that of the avail- 
able SDSS images. It was given by the Poisson noise associated 
to the photon counts due to both the galaxy and sky background 
and read-out noise (RON) of the CCD. The pixel scale, CCD 
gain and RON were 073946 arcsec pixel -1 , 4.72 e~ ADIT 1 , and 
5.52 e~, respectively. They mimic the instrumental setup of the 
SDSS images. In order to account for seeing effects, the images 
of the artificial galaxies were convolved with a Moffat PSF with 
FWHM=2.77 pixels and /3 = 3.05 (see Sect.©. 

The artificial galaxies do not match the sample galaxies since 
we did not account for their redshift distribution. Due to this 
issue, the fractions given in Tab. [3] and |4] do not represent at 
all estimations of the absolute bar fraction lost in real galaxies. 
However, they are useful to test the efficiency of the two pro- 
posed methods for detecting bars, in order to fine tune their free 
parameters and understand possible biases in the results. The ar- 
tificial galaxies were sized in pixels. In this way, the images are 
somewhat 'dimensionless' and the performances of two meth- 
ods can be assessed by converting the scale lengths from pixel to 
physical units according to the distance of the objects. In order 
to cover the full range of bar lengths and redshifts, the shortest 
bars (r bar = 0.5 kpc) were scaled assuming a redshift z = 0.04, 
while the largest ones (r bar = 5 kpc) were placed at z = 0.01. 

4.2. Testing the ellipse fitting method 

The ellipse fittin g meth od has two free parameters, Ae and 
APA. lLaine et all d2002l) a dopted Ae = 0.1 and APA = 20° , 
IMarinova & Jogedl20a7l) , iMenendez-Delmestre et alj d2007l) . 
and lBarazza et all d2008l) adopted Ae = 0.1 and APA = 10°. 

We applied the ellipse fitting method to the artificial galaxies 
by adopting Ae = 0.1, 0.08, 0.05 and APA = 10°, 20°, 30° to find 
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Table 3. Percentage of galaxies classified as barred and unbarred 
galaxies erroneously found to be barred by applying the ellipse 
fitting method to the sample of artificial galaxies. 



APA 




Barred 




Unbarred 




Ferrers Bars 


Freeman Bars 


Flat Bars 




Ae = 0.10 


10° 


39 


46 


22 


2 


20° 


53 


58 


34 


3 


30 s 


57 


63 


37 


4 


Ae = 0.08 


10° 


48 


55 


30 


3 


20° 


62 


67 


40 


5 


30° 


66 


71 


43 


6 


Ae = 0.05 


10° 


65 


70 


43 


9 


20° 


77 


81 


54 


13 


30° 


79 


83 


57 


15 



the best combination of Ae and APA maximising the bar identi- 
fications and minimising the bad and/or spurious detections. The 
results are given in Table[3] 

The flat and Freeman bars are the most difficult and eas- 
iest bars to be detected, respectively. This means that the el- 
lipse fitting method detects more efficiently the bars with sharp 
ends than those characterised by a smooth transition to the disc. 
Moreover, adopting APA = 20° instead of APA =10° increases 
the fraction of bar detections by 10% for all the bar types, while 
the increment between APA = 20° and APA = 30° is only about 
4%. 

Spurious detections correspond to unbarred galaxies which 
are erroneously found to be barred. In order to estimate their 
fraction, we applied the method to the sample of unbarred ar- 
tificial galaxies we built to this aim. The results for the differ- 
ent values of Ae and APA are also given in Table [3] Bad detec- 
tions correspond to barred galaxies for which we obtained a bad 
measurement of the bar length. In order to estimate this fraction 
(which is not reported in Table [3]), we compared the bar lengths 
known for the artificial galaxies with the rf derived by applying 
the ellipse fitting method. We derived the median and standard 
deviation of the relative error between the known and measured 
bar length using a 3<x clipping iterative procedure. We consid- 
ered as bad detections the measurements with a relative error 
larger than 3<x with respect to the median. 

The fraction of galaxies classified as barred versus the 
bad/spurious detections are shown in Fig. [3] as a function of Ae 
and APA. The optimal configuration is Ae = 0.08 and APA = 20° 
since the fraction of detections increases by more than 10% with 
respect to Ae = 0.1 and APA = 10°, while the fraction of 
bad/spurious detections is always lower than 10%. It is worth 
noticing that for A e = 0.08 and APA = 30° the fraction of 
bar detections rises by about 3% with respect to A e = 0.08 and 
APA = 20°. But the fraction of bad/spurious detections increases 
too. For example, for the Ferrers bars such a fraction is even 
larger than 10%. 

The bar lengths we measured as and rPj* in the artifi- 
cial galaxies are shown in Fig. [4] The bar length is underesti- 
mated when rf is used, as found bv lMichel-Dansac~& Wozniak 
(2006) too. This is particularly true for the Ferrers bars where the 
measured bar lengths are 51% shorter than the real ones. In con- 
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Fig. 3. Fraction of artificial galaxies correctly classified as barred 
vs fraction of bad/spurious bar detections using the ellipse fit- 
ting method with APA=10° (asterisks), 20° (filled circles), and 
30° (squares) and Ae =0.1 (black symbols), 0.08 (red symbols), 
and 0.05 (blue symbols). The large, medium, and small symbols 
correspond to Ferrers, Freeman, and flat bars, respectively. 

trast, Freeman and flat bars were better determined, their mea- 
surements being shorter by 30% and 19%, respectively. The bar 
length is underestimated when r?^ is used for the Ferrers bars 
(11%), but it is overestimated for the Freeman (8%) and flat bars 
(28%). These results show the possibility of define an empirical 
correction to the bar length, knowing the bar type in advance. 

4.3. Testing the Fourier analysis method 

The deprojection of the galaxy image is a crucial step in applying 
the Fourier method, which has two free parameters, A(/2//o) and 
A(f>2- The inclination and major-axis position angle of the galaxy 
disc can be obtained b y either fitting ellip se to the outermost 
galaxy isophotes (e.g., lAguerri et all 12003b or minimising the 
m = 2 Fourier m ode in the outermost regions of the galaxy (e.g., 
iGrosboll |T985). We applied these two methods to the artificial 
galaxies and found that ellipse fitting gave lower errors (about 
3°) on both /and PA. 

We applied the Fourier method to the artificial galaxies by 
adopting A(/ 2 // ) = 0.2, 0.1, 0.08, 0.05 and A<p 2 = 10°,20°,30° 
to find the best combination of A(/2//o) and Afc, maximising 
the bar identification and minimising the bad and/or spurious 
detections. The results are given in Tableland in Fig. [5] 

In general, the Fourier method is less efficient in detecting 
bars than the ellipse fitting method. We found that A(l2/Io) = 0.2 
has to be adopted to have a fraction of bad/spurious detections 
lower than 10%. We also adopted A<f>2 = 20° because it increases 
detections by more than 10% and gives less bad/spurious detec- 
tions with respect to A02 = 30°. The bar lengths we measured as 
r bai Ul er ' n t ' le artificial galaxies are shown in Fig. [4] The method 
recovers the bar length with the best accuracy for the Ferrers 
bars (3% error), while the bar length is over estimated for the 
Freeman (28%) and flat bars (46%). 

5. Bar fraction 

Both the ellipse fitting and Fourier method were applied to our 
sample of 2106 disc galaxies. We found that the fraction of 
galaxies classified as barred depends strongly on the technique 
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Fig. 4. Bar length measured with the maximum ellipticity (top panels), position angle (middle panels), and Fourier analysis method 
(bottom panels) for Ferrers (left panels), Freeman (central panels), and flat bars (right panels) in artificial galaxies. Mean relative 
error (defined as the difference between input and measured values) and standard deviation for the measurements are given in each 
panel. 



Table 4. Percentage of galaxies classified as barred and unbarred 
galaxies erroneously found to be barred by applying the Fourier 
analysis method to the sample of artificial galaxies. 



A0 2 




Barred 






Unbarred 




Ferrers Bars 


Freeman Bars 


Flat Bars 








A(/ 2 //o) = 


0.20 






10° 


29 


26 




18 
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20° 


44 


42 




31 
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30° 


51 


51 




37 
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A(/ 2 //o) = 


0.10 






10° 
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35 
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20° 
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52 
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30° 


69 
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31 
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40 


22 
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67 




58 
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0.05 
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53 


29 


20° 


75 
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70 


41 


30° 


81 
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48 



adopted for the analysis: it is 45% with the ellipse fitting method 
and 26% with the Fourier method. Although the Fourier method 
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Fig. 5. Fraction of artificial galaxies correctly classified as barred 
and fraction of bad/spurious bar detections using the Fourier 
analysis method with A02 = 10° (asterisks), 20° (filled cir- 
cles), and 30° (squares) and A(/2//o) = 0.2 (green symbols), 0.1 
(black symbols), 0.08 (red symbols), and 0.05 (blue symbols). 
The large, medium, and small symbols correspond to Ferrers, 
Freeman, and flat bars, respectively. 



was demonstrated to be less efficient than ellipse fitting in de- 
tecting bars, this difference is larger than that expected from the 
analysis of the artificial galaxies. To investigate this issue, we 
took into account the morphological classification of the galax- 
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Table 5. Median values of the bar radius for the different galaxy 
types. 
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Fig. 6. The r-band images (left panels), ellipticity (central pan- 
els), and m — 2 Fourier amplitude radial profiles (right panels) 
of the early-type spiral SDSSJ03 1947.01 +003504.4 (top panels) 
and the late-type spiral SDSSJ020159.33-081441.9 are anal- 
ysed by both methods (bottom panels). For each galaxy the ver- 
tical dotted line corresponds to the value of . 



ies found to be barred. According to the ellipse fit method 29%, 
55%, and 54% of the lenticular, early-type and late-type spiral 
galaxies, respectively, are barred. They are 29%, 33%, and 17%, 
respectively, with the Fourier method. Therefore, both methods 
obtained a similar of fraction of barred lenticular galaxies, while 
the Fourier method is less efficient in detecting bars in spirals, 
and particularly in late-type spiral galaxies. 

An example is shown in Fig. [6] The early-type spi- 
ral SDSSJ03 1947.01 +003504.4 and the late-type spiral 
SDSSJ020159.33-081441.9 are analysed by both methods. The 
bar of the early-type spiral was detected by both methods and 
the measured bar lengths are in agreement. In fact, the radial 
profiles of e and h/Io show a local maximum at about 5", 
where the PA and <p2 are constant. On the contrary, the bar of 
the late-type spiral was detected only by the ellipse fit method. 
The radial profiles of e and hlh show a local maximum at 
different radii (about 10")- The hlh maximum is located in 
the spiral arm region, where <p2 is not constant. Therefore, the 
bar of this galaxy was not detected by the Fourier method. We 
conclude that bars with sharp ends are detected by both ellipse 
fitting and Fourier methods. But, the bars of galaxies with lenses 
or strong spiral arms are more easily detected with the ellipse 
fitting method. This kind of bars is usually found in late-type 
spirals. These large differences in the bar fractions between the 
two methods could bias our conclusions. For this reason, we 
will study the photometrical parameters of the bars by adopting 
only the ellipse fitting method. 

As an additional check, we visually classified all galaxies in 
barred and unbarred. The visual classification was done by two 
of us (JALA and JMA). Both the classifications were in close 
agreement and only their mean is reported. No attempt to clas- 
sify the galaxies according to their Hubble type was done. The 
difference between the bar fractions found with the visual and 
the automatic classification (~ 7%) is our best estimate of the 
fraction of undetected bars in the galaxy sample. We obtained 
that the global fraction of barred galaxies in our sample was 
38%. Taking into account the different morphological types we 
found that 22%, 52% and 48% of the lenticular, early-type and 
late-type spiral galaxies, respectively, were barred. 
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Our bar fraction (45%) of disc galaxies in the local uni- 
verse is in goo d agreement with recent r esults obtained in op- 
tical bands by IMarinova & Jogeel d2007l 44% in the B band), 
iReese et alj d2007l 47% in the / band). and lBarazza et alj (2008, 
50% in the r band). In addition, our finding that early- and late- 
type spirals host a larger fracti on of bars than lentic ular galax- 
ies, was also in agreement with Barazz a"et alj ((2008) since they 
found that disc-dominated galaxies with low bulge-to-disc lumi- 
nosity ratio display a hig her bar fraction than galax ies with sig- 
nificant bulges. However. IMarinova & Jogeel d2007l) did not find 
any difference in the bar fraction in the NIR as a function of the 
Hubble type, we argue that their result is biased by their smaller 
coverage of the Hubble sequence since neither lenticulars nor 
Sd/Sm galaxies were taken into account. The same considera- 
tion can be applied to the results found by Knapen et al. (2000) 
and lEskridge et al.1 d2000t) . 



6. Bar properties 

6.1. Bar length 

The distributions of the bar lengths and normalised bar lengths 
were derived for the sample galaxies after deprojection on the 
galaxy plane. They are shown in Fig. [7] Both rf and rP^ were 
considered, and the galaxy size was defined as r ga i = 2 x r p , 
where r p is the Petrosian radius from SDSS. The median values 
derived for the different morphological bins are given in Tab. [5] 
The values of rf are systematically smaller than those of rj^., as 
expected from the measurements of the artificial galaxies. The 
comparison of our results, with previous works where the rf^ 
values are reported, gives us a good agreement. For example, 
Erwin (2005) found a median bar length of 3.3 kpc, Marinova 
& Jogee (2007) calculate a mean value of 4 kpc and Menendez- 
Delmestre et al. (2007) obtain a median value of 3.5 kpc. These 
results hold even considering only the bars with a length larger 
than 2 kpc. This limit corresponds to minimum bar length we are 
able to resolve all throughout our range of distances. 

Since the bar length is strongly dependent on the method 
used to derive it, we can not conclude much about the correlation 
between the bar length and the morphological type. According to 
rf , the lenticulars host the shortest bars, while according to r? A 

bar' ' & bar 

their bars are the longest ones. As far as the median bar length of 
the spirals concerns, the late-type spirals host shorter bars with 
respect to the early-type ones (Tab. |5). Nevertheless, a correla- 
tion between the bar length and galaxy size was found. Thus, 
larger bars are located in bigger galaxies (Fig. [8]). The correla- 
tion is independent of the adopted method to measure the bar 
length. It holds for the different morphological bins too, being 
stronger for late-type spirals (r = 0.52) and weaker for SO galax- 
ies (r = 0.38). This relation could indicate a link between the 
formation and evolution processes between of bars and galaxy 
discs. A similar correlation was found by IMarinova & J ogee 
( 2007), although a quantitative comparison with them is not pos- 
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Fig. 7. Distribution of the bar length (left panels) and normalised bar length (right panels) in lenticulars (full line), early-type spirals 
(dotted line), and late-type spirals (dashed line). The bar length was measured with both the maximum ellipticity (top panels) and 
PA method (bottom panels). Arrows mark the median values of the distributions. 
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Fig. 8. Galaxy radius r ga i versus bar length measured with the 
PA method (top panel) and maximum ellipticity method (bottom 
panel). The solid line represents the linear regression through all 
the data points. The Pearson correlation coefficient (r) and the 
result of the linear fit are given. 



sible due to the different band-passes and different definition of 
the galaxy radius they adopted. 



6.2. Bar strength 

The bar strength represents the contribution of the bar to 
the total galaxy potential. Several methods have been devel- 



oped to measure it (see [Laurikain en et all 120071 and references 
therein). Nowadays, the most commo nly used parameter m ea- 
suring the bar strength is Q % defined bv lButa & Blockl (120011) . It 
can be accu rately estimated by analysing near-infrared image s 
dButa et all [2003t iBlock eTail l2004t iLaurikainen et all 120071) . 
which are not av ailable for our sample galaxies drawn from the 
SDSS. However, lAbraham & Merrifieldl (120001) defined another 
bar strength parameter given by 



/bar = - (arctan(l - e bar )" 1/2 - arctan(l - e bal -) +1/2 ) , 



(29) 



where eb ar is the bar ellipticity me asured at rf ■ It correlates with 
Q g (e.g., Lauri kainen et alll2007l) and was adopted for our sam- 
ple galaxies. 

We did not adopted any minimum value for the bar ellipticity. 
We found e ba i,min = 0.16, which is close the minimum ellipticity 
adopt ed in other studies (e.g., ebar.min = 0.2. lMarinova & JogeeL 
2007). The distributions of the bar strengths we derived for the 
different morphological types are shown in Fig. [9] The median 
values for the bar strengths of the lenticular, early-type, and 
late-type spiral galaxies are 0.16, 0.19, and 0.20, respectively. 
Indeed, we found a significant difference between the lenticu- 
lar and spiral galaxies. They are characterised by different dis- 
tributions, as confirmed at a high confidence level (> 95%) by 
a Kolmogorov-Smirnov ( KS) test. Using four differ ent methods 
to derive the bar strength, Laurikaine n et all (120071) also found 
that SO galaxies host significantly weaker bars than the re st of 
disc galaxies, t his re sult was hold also by dDas et all 120031) and 
Barazzaet al. (2008) using only the bar ellipticity. In contrast, 
Mari nova & Joged d2007l) found that the ellipticity of the bar is 
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Fig. 9. Distribution of the bar strength in lenticulars (full line), 
early-type spirals (dotted line), and late-type spirals (dashed 
line) for the whole sample (top panel) and for the galaxies with 
'"bar > '"bd (bottom panel). Arrows mark the median values of the 
distributions. 



Fig. 10. Cumulative distribution of barred (solid black line) and 
unbarred (dashed black line) galaxies as function of the local 
galaxy density. Cumulative distribution of barred (long-dashed 
blue line) and unbarred (dashed-dotted blue line) galaxies af- 
ter excluding those at less than 7 Mpc from the nearest edge of 
SDSS is also overplotted (see text for more details). 



practically independent of the Hubble type. But they not con- 
sider SO galaxies which are those actually making the difference. 

However, it could be possible that the presence of a large 
bulge could affect the measurement of the bar ellipticity, and 
therefore the calculation of the bar strength. In order to address 
this issue, we performed a further test. We fitted an exponential 
law to the outer parts of the surface-brightness profiles of our 
barred galaxies. Then, we computed the radius at which the 
galaxy surface brightness profile exceeds the fitted exponential. 
This radius represents an estimate of the extension of the region 
where the bulge contribution dominates the light of the galaxy. 
At this point, we selected a subsample of barred galaxies with 
1>ar > /"bd- We recalculated the mean strength of these bars by 
splitting the sample in lenticulars, early- and late-type spirals. 
As expected, in the new subsample of galaxies we lost the weak- 
est bars, especially in the lenticular galaxies. Nevertheless, the 
final result is the same: the lenticulars have weaker bars than the 
early- or the late-type spirals. 

Some numerical simulations of bar formation and evolution 
propose that bars can be formed and destroyed fastly due to 
the accretion of gas towa rds the central regions of the galax- 
ies dPfenniger & Normanl |1990t iBournaud & Combesl 120021; 
Bournaud et all 120051) . In this framework, due to the fast de- 
struction and re-formation of bars, we would expect a bimodal 
distribution of the bar strength at least for gas-rich galaxies like 
sthe late-type barred ones. The absent of this bimodality in the 
bar strength for all galaxy types showed in Fig. [9] could be 
against those bar formation and evolution scenarios. 



7. Bars and galaxy properties 

7.7. Galaxy local environment 

Due to our selection criteria we excluded all the strongly dis- 
turbed/interacting galaxies. Nevertheless, we calculated for each 
sample galaxy the loca l density following the prescriptions of 
Balogh et al. (2004a b) in order to investigate the relation be- 
tween the bar properties and local environment of the host 
galaxy. The number density of local galaxies was computed us- 



ing the distance of the galaxy to its fifth nearest neighbour 
galaxy. Thus, a projected galaxy density could be defined as 

nd 1 ^ 

This was computed with those galaxies located in a velocity 
range of +1000 km s _1 from the target galaxy to avoid back- 
ground/foreground contamination. For sample galaxies with- 
out a measured redshift we imposed a luminosity constraint of 
+3 mag to derive the galaxy density as done by Balog fTet al.l 
(l2004aibl) . 

Fig. [TUl shows the fraction of barred and unbarred galaxies in 
our sample as function of the local galaxy density. In the range 
of galaxy density covered by our sample, there is no relation 
between the presence of a bar and the environment of the host 
galaxy. The same is true even if the galaxies of different mor- 
phological type are considered independently. In addition, we 
did not find any correlation between the bar length or strength 
and the local galaxy density. 

However, a series of caveats must be taken into account when 
the distance neighbour method is applied. 

The limited area of sky covered by the survey implies that 
usually the estimated density is lower than the true one. In fact, 
due to edges and holes in the survey, the value of d 5 could 
be overestimated and the derived density underestimated. Two 
methods were applied to our sample in order to test the robust- 
ness of our result against the edge effects. The first consists in 
removing all the galaxies whose distance to the su rvey edges is 
smaller than the measured fi fth neighbour distance dMiller et al.L 
l200l iBalogh et all l2004bb . In this way, we ensure that the re- 
maining galaxies have an unbiased estimation of the local den- 
sity. The fraction of galaxies we measured in each density bin 
after applying such a correction is shown in Fig.Qj] It is worth 
noticing that no galaxy is available in the bin of lowest local den- 
sity (S5 < 0.01 Mpc -2 ). This means that such a low density can 
not be derived for any galaxy of the sample. On the contrary, the 
local density is correctly derived for all the sample galaxies with 
E 5 > 1 Mpc -2 . 
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Unfortunately, this method biases the distribution toward 
over-dense environments. This bias can be reduced by exclud- 
ing all the galaxies within a given distance to the nearest edge 
(Coo per et al.ll2005l) . We excluded the galaxies at less than 2, 5, 
and 7 Mpc from the nearest edge of SDSS (Fig.lTO. If the lowest 
density bin is excluded, the method introduces only a weak con- 
tamination toward high density environments. An optimal given 
distance of 7 Mpc was found by calculating the maximum dis- 
tance to the fifth neighbour in the bin of lowest density. In this 
way, the local density measurements are reliable and there is no 
bias toward the high-density environments, i.e. all the bins have 
almost the same number of galaxies. Fig. QJJshows the fraction 
of barred and unbarred galaxies in our sample as function of the 
local galaxy density after excluding galaxies at less than 7 Mpc 
from the nearest edge of SDSS. They are about 40% of the total. 
We do not find any difference between the environment of barred 
and unbarred galaxies. 

Also the selected redshift range could lead to underestimate 
the local density. We circumvent this problem by defining a new 
volume-limited sample in a wider redshift range (0 < z < 0.06) 
taking into account the adopted velocity range. 

Finally, the local density distribution could be biased by 
SDSS fiber collision which prevents to measured galaxies closer 
than 55" with respect to each other. In SDSS-DR5 the net effect 
of fiber collision is a loss of 6% of the galaxies in the photometric 
catalogue that would othe rwise be in the spectroscopic catalogue 
dCowan&Ivezicl [2008). In our case, this value represents an 
upper limit since all the galaxies in photometric catalogue were 
taken into account in calculating the local density. 

According to numerical simulations, galaxy mergers and 
interactions are mechanism s which should drive the for- 
mation of bars dGerin et all Il990t iMiwa & Noguchl Il998t 
iMastropietro et all 20051) . Therefore, we could expect that frac- 



tion of barred galaxies increases with the local density. But, the 
observational proofs about the influence of the envir onment on 
bar fo rmation and evolution are few. For example, Thompson 
(1981) suggested a link between bar formation and local galaxy 
environment by observing that the fraction of barred galaxies in- 
creases towards the core of Coma cluster. But, this is not case 
in the wide range of densities we explored (0.01 < E5 < 100 
Mpc 2 , Fig. [Toll. For the lowest density bin o ur fraction of barred 
galaxies is even smaller than 60% found bv lVerlev et ail (120071) 
by analysing the optical images of 45 isolated galaxies. Recently, 
( Marino v a et al. L 120 08) shows that the cluster environment does 
not strongly affect to the bar fraction. We argue that for non- 
interacting and undisturbed galaxies the environment do not play 
a major role in the formation and evolution of their bar. 

Fig. [TO] shows that 80% of the sample galaxies are located 
in very low-density environments (S5 < 1 Mpc 2 ). The lo- 
cal density of the remaining 20% (corresponding to more than 
400 galaxies) covers mostly typical values measured for loose 
(£5 > 1 Mpc -2 ) and compact galaxy groups (S5 ~ 10 Mpc -2 ). 
Nevertheless, the fraction of barred galaxies does not depend 
on the local density also for these galaxies. Therefore, we con- 
clude that the environment does not play an important role in 
the formation of bars, at least over the observed range of lo- 
cal densities. Moreover, it does not account for the variation of 
the central light concentration and galaxy colours discussed in 
Sects. [7721 and 17.31 respectively. Similarly, low density environ- 
ments, as those reported here, do not also account for variations 
in other galaxy prop erties, such as the blue galaxy fraction (e.g., 
Agu erri et al.L 120071) . However, we can not infer that close in- 
teractions do not affect bar formation and evolution, because we 
selected only non-strongly disturbed/interacting galaxies. 
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Fig. 11. Fraction of galaxies remaining in the sample after cor- 
recting for survey edge effects by removing all the galaxies with 
an edge of the survey closer than the measured fifth neighbour 
distance (blue squares), 2 Mpc (black diamonds), 5 Mpc (green 
triangles), and 7 Mpc (red asterisks). Poissonian errors are given. 
Horizontal lines show the width of the local density bins. 



This result i s in ag r eement with the numerical simulations 
by iHeller et all (12007b : iRomano-Diaz et al.1 (120081) . who argue 
that there is no difference between the bar fraction for field and 
cluster galaxies. They claim that the bar evolution is mainly 
driven by the dark matter subhalos, which surround all the bright 
galaxies and do not depend of their environment. These subhalos 
could host faint galaxies, which are not visible in our images. 



7.2. Central light concentration 

Fig. [12] shows the distribution of the C parameter for the barred 
and unbarred sample galaxies and their cumulative distribution 
functions. Both types of galaxies are characterised by different 
distributions as confirmed at a high confidence level (> 99% ) 
by a KS test. It is worth noting that the differences between both 
distributions is due to galaxies with higher central light concen- 
tration. This result holds even if we take into account the con- 
tamination of ellipticals into our sample of disc galaxies (Sect. 

0. 

iBarazza et al. I (120081) found that the fraction of barred galax- 
ies is higher for the galaxies with a smaller value of the Seisic 
parameter (i.e., the less-concentrated galaxies). We confirm their 
findings. Fig. [12] shows that the number of barred and unbarred 
galaxies is clearly different for galaxies with high values of C, 
being the fraction of barred galaxies smaller than the unbarred 
ones. 

Since light concentration is correlated with the cen- 
tral velocity dispersion, the previous result implies that, in 
some way, galaxies with higher central mass concentrations 
tend to inhibit the formation and/or evolution of bars. This 
is in agr eement with the results of the numerical experi - 
ments by iPfenniger & Normanl d 19901): iNorman et al.l dl 996 ) ; 
lAthanassoula & Misiriotisl d2002l) : iDeb attista et al. (2006), who 
showed that the presence of a large bulge weakens the bar. 
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Fig. 12. Number of barred (full line) and unbarred disc galaxies 
(dotted line) as function of the light concentration (upper panel). 
Cumulative distribution of barred (full line) and unbarred disc 
galaxies (dashed-line) as function of the light galaxy concentra- 
tion (lower panel). 




Fig. 13. Cumulative distribution of barred (solid line) and un- 
barred (dashed line) as function of the g — r galaxy colour. 



7.3. Galaxy colours 

Fig. [T3l shows the cumulative distribution functions of the g — r 
colour for the sample galaxies. Barred and unbarred galaxies are 
characterised by different distributions as confirmed at a high 
confidence level (> 99%) by a KS test. Thus, barred galaxies 
are bluer than unbarred ones. We can explained this effect as due 
to the larger fraction of barred galaxies observed in the late-type 
systems, which are systematically bluer than the early-type ones. 
Similar colour differ ence between barred and unbarred galaxies 
was also reported bv lBarazza et alj (120081) . 



8. Conclusions 

We have studied the fraction and properties of bars in a sam- 
ple of 2106 disc galaxies extracted from the SDSS-DR5. This 
is a volume-limited sample of undisturbed and non-interacting 
galaxies in the redshift range 0.01 < z < 0.04, with an absolute 
magnitude M r < -20, and an inclination i < 60°. 



The sample galaxies have been classified lenticulars, early- 
(i.e, SO/a - Sb) or late-type spirals (i.e, Sbc - Sm) according 
to their light concentration. The light concentration was defined 
as C — rgo/rgo, where r$o and ryo are the radii enclosing 50% 
and 90% of the total galaxy light, respectively. These radii are 
available in the SDSS database for all objects of our sample. 
The correlation between the light concentration and morpholog- 
ical type was derived from a subsample of 612 galaxies listed 
in RC3, whose morphological classification was already known. 
We found that the fraction of lenticulars, early- and late-type 
spirals among the selected disc galaxies is 39%, 28%, and 33% 
respectively. 

We derived the fraction of barred galaxies by analysing the 
SDSS r-band images with the ellipse fitting and Fourier analysis 
methods. They consist in looking for a local maximum in the ra- 
dial p rofile of ellipticity ( associated to a constant position angle; 
e.g.. IWozniak et all 1 19951) and m - 2 relative F ourier compo- 
nent (a ssociated to a constant phase angle; e.g., lAguerri et all 
2000a), respectively. 

The bar fraction depends strongly on the technique adopted 
for the analysis. By extensive testing on a large set of artificial 
galaxies, we concluded that the both methods are efficient in de- 
tecting the bars with sharp ends, such as the Ferrers and Freeman 
bars. On the contrary, the flat bars, which are characterised by a 
smooth transition to the disc, are more easily detected by the 
ellipse fitting method. The ellipse fitting method is more effi- 
cient in detecting bars in galaxies with lenses and spiral arms, 
where the m = 2 relative Fourier component shows multiple 
maxima (but it is not associated to a constant phase angle). This 
is the case for late-type spiral galaxies. For this reason, we de- 
cided to rely only onto the results obtained with the ellipse fitting 
method. We found that 45% of the selected disc galaxies host a 
bar in agreement with previous findings in o ptical wavebands 
dMarinova & Joged 120071; iReese et all 120071) . The fraction of 
bars in the three different morphological bins is 29%, 55%, and 
54% for lenticulars, early- and late-type spirals, respectively. By 
classifying visually the galaxies in barred and unbarred we ob- 
tained similar bar fractions as those reported by the ellipse fitting 
method. 

The bar length was obtained by measuring the radius r? at 
which the maximum ellipticity was reached and as the radius 
r? A at which the PA changes by 5° with resp ect to the value 
corre sponding to the maximum ellipticity (e.g.. IWozniak et all 
Il995l) . According to the analysis of the artificial galaxies, the 
bar length is underestimat e d wh en rf is used, as found by 
Mich el-Dansac & Wozniakl ([2006) too. This is particularly true 
for the Ferrers bars where the measured bar lengths are 51% 
shorter than the real ones. In contrast, Freeman and flat bars were 
better determined, their measurements being shorter by 30% and 
19%, respectively. The bar length is underestimated when is 
used for the Ferrers bars (11%), but it is overestimated for the 
Freeman (8%) and flat bars (28%). These results show the possi- 
bility of define an empirical correction to the bar length, knowing 
the bar type in advance. We obtained that the bar lengths (when 
normalised by the galaxy size) are larger in lenticular galaxies 
than those presented in early- and late-type ones. This finding 
is independent of the method used for measuring the bar length, 
and statistically significant according with the KS test. We also 
found a correlation between the bar length and galaxy size. This 
correlation is also independent of the method used for measur- 
ing the bar length. It holds for the different morphological bins, 
being stronger for late-type spirals (r = 0.52) and weaker for 
lenticular galaxies (r = 0.38). The larger bars are located in 
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larger galaxies, indicating an interplay between the bar and disc 
in galaxy evolution. 

The b ar strength /bar was estimated following the parametri- 
sation by Abra ham & Merrifieldl yOOO) which requires the mea- 
surement of the bar ellipticity. The median values for the bar 
strengths of the lenticular, early-type, and late-type spiral galax- 
ies are 0. 16, 0. 19, and 0.20, respectively. The bars of the lenticu- 
lar galaxi es were found to be wea ker than those in spirals, as 
found by lLaurikainen et al.l d2007l) too. The fact that the bar 
strength distribution are unimodal for all galaxy types argues 
against evolutionary models in which bars would be formed and 
destroyed in short timescales. 

No difference between the local galaxy density was found 
between barred and unbarred galaxies in our sample. Thus, the 
local environment does not seem to influence bar formation. 
Moreover, neither the length nor strength of the bars are cor- 
related with the local galaxy environment. The previous results 
are even true for the subsample of our galaxies located in the 
more dense environments (logl^) > Gal/Mpc 2 . Those envi- 
ronments could be similar to those showed by galaxy groups or 
weak clusters of galaxies. Thus, we can say that even for the 
densest environments, the global environment do not play an im- 
portant role in the bar formation. However, we can not exclude 
than close galaxy-galaxy encounters would trigger the bar for- 
mation, as they were excluded from our sample. These results 
indicate that formation and evolution of the bars in the studied 
sample depend mostly on internal galaxy processes rather than 
external ones. 

A statistical significant difference between the central light 
concentration of barred and unbarred galaxies was found. The 
bars are mostly located in less concentrated galaxies. This 
difference could explain the lower fraction of bars detected 
in SO galaxies with respect to spirals. Since the SO galax- 
ies host weaker bars than spirals, we conclude that central 
light concentration is an important factor driving the bar for- 
mation an d evolution. In fact, accord i ng to the numerica l 
simulations iPfenniger & Nqrmanl (ll990h:lNorman et al.l (Tl 996): 
lAthanassoula & Misiriotisl d2002l) ; iDebatfista et all (120061) . the 
bars are weakened by large bulges. Finally, bars are mainly 
hosted by bluer late-type spirals. We argue that this is due to 
late-type galaxies have larger bar fraction than early-type ones. 
Similar results were also found in previous works as Barazza et 
al. (2008). 

The sample of galaxies presented in this study is one of the 
largest samples presented in the literature. This large number 
of studied galaxies makes that the conclusions reported in the 
present work about the different observational bar properties is 
stronger than those obtained with smaller number of galaxies. 
For this reason, the present work will be useful for constraining 
future theoretical works about formation and evolution of bars 
in disc galaxies. 
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